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Abstract 

Background: Multi-nucleon transfer reactions have recently attracted attention as a possible 
path to the synthesis of new neutron-rich heavy nuclei. 

Purpose: We study transfer reactions involving massive nuclei with the intention of under¬ 
standing if the semi-classical model GRAZING coupled to an evaporation and fission competition 
model can satisfactory reproduce experimental data on transfer reactions in which fission plays a 
role. 

Methods: We have taken the computer code GRAZING and have added fission competition to 
it (GRAZING-F) using our current understanding of F^/F f, fission barriers and level densities. 

Results: The code GRAZING-F seems to satisfactory reproduce experimental data for +lp, 
+2p and +3p transfers, but has limitations in reproducing measurements of larger above-target 
and below-target transfers. Nonetheless, we use GRAZING-F to estimate production rates of 
neutron-rich N = 126 nuclei, actinides and transactinides. 

Conclusions: The GRAZING code, with appropriate modifications to account for fission decay 
as well as neutron emission by excited primary fragments, does not predict large cross sections for 
multi-nucleon transfer reactions leading to neutron-rich transactinide nuclei, but predicts oppor¬ 
tunities to produce new neutron-rich actinide isotopes. 

PACS numbers: 25.70.Hi,25.60.Je,25.85.-w,24.10.Lx 
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I. INTRODUCTION 


Experimentalists have had a long-standing interest in multi-nucleon transfer reactions [U 
[2] hoping to synthesize new neutron-rich isotopes not normally accessible by neutron capture 
and fusion reactions [3H9]. Cross sections of actinides produced in transfer reactions using 
light and heavy projectiles and actinide targets were measured by the chemical separation 
of the products in a series of experiments in the late 70’s and 80’s. The systematic trend 
that emerged after the series of experiments with U, Cm and Cf targets is that the use of 
transfer reactions to produce unknown neutron-rich actinides is favorable for below-target 
species and limited for above-target species. The production of neutron-rich trans-target 
nuclides up to Fm and Md with cross sections ~ 0.1/ib were observed. The basic problem 
in making heavier nuclei was that the higher excitation energies that led to broader isotopic 
distributions caused the highly excited nuclei to fission. 

The interest in transfer reactions has been recently boosted by the prediction of larger 
than expected cross sections for the production of heavy nuclei, within the framework of a 
dynamical model based on Langevin equations, by taking advantage of shell effects which 
may favor a large flow of nucleons resulting in the formation of surviving heavy nuclei 
pm ITT] . In this picture, low-energy multi-nucleon transfer reactions of very heavy nuclei, 
such as U-|-Cm, may produce one primary reaction product in the vicinity of Z = 82, 
N = 126 closed shells, leaving the second primary product in the actinide or transactinide 
region with very low excitation energy and thus, with increased probability of surviving 
hssion. This model was able to account for the previously measured radiochemical data 

ra- 

The motivation and interest in multi-nucleon transfer reactions in Ref. [lO] and the 
present paper is two-fold: (a) the possibility of producing the most neutron-rich heavy 
nuclei for studies using nuclear spectroscopy, atomic physics and chemistry and (b) the 
difficulty in pursuing the study of nuclei with high atomic numbers using fusion reactions. 
Traditional “cold” fusion reactions have production cross sections of 10 — 100 fb beyond 
Z = 112, and “hot” fusion reactions have cross sections of the order of a few pb for elements 
Z ~ 118. The upper limit cross sections for Z = 119 and Z = 120 have been established to 
be of the order of 100 fb |T3]. This difficulty has spurred the renewed interest in low-energy 
multi-nucleon transfer reactions as a way of accessing new neutron-rich transactinide nuclei 
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that are closer to the “island” of stability near the nentron shell N = 184 not accessible by 
fnsion reactions. 

Mnlti-nncleon transfer reactions in the qnasi-elastic and deep-inelastic regimes have been 
extensively modeled with the semi-classical description of Winther [lH [15], implemented in 
the compnter code GRAZING [T6|. This code considers the mnlti-step exchange of nncleons 
between the colliding nnclei in classical trajectories calcnlated with a Gonlomb plus nuclear 
interaction. GRAZING is known to have shortcomings, i.e. the initial deformations of the 
nuclei is not taken into account and neutron evaporation from the primary products is the 
only de-excitation mode considered. As a result, the code has mainly been used to predict 
yields in light projectile reactions with medium to heavy targets in which the hssility of 
the reaction products studied is small [IM5|. The theoretical formalism of GRAZING is 
described in depth by Winther [TTl |15]. An outline of the main ingredients and approxi¬ 
mations of the model can be found in the topical review by Gorradi, Pollarolo and Szilner 
[2n|- Multi-nucleon transfer reactions have also been studied theoretically using the Fokker- 
Plank equation [27], the hnite-range DWBA model |2B], the Di-nuclear System model [22] . 
the time-dependent Hartree-Fock theory ED], and the Langevin equations [TO] . 

The GRAZING code has recently been informally used to predict yields of products in 
reactions with planned radioactive beams (EURISOL) and isotope “factories” (GARIBU), 
in some cases with actinide targets. In this paper we present an extension to GRAZING in 
which not only neutron evaporation from the excited primary products is considered, but 
also hssion competition. With such additions to the code, reactions where hssion effectively 
competes with neutron emission, e.g. the U-I-Gm reaction, can be studied and compared to 
experimental data and other models. 


II. NEUTRON EVAPORATION AND FISSION COMPETITION 


The competition between neutron emission and hssion is simulated with the classical 
formalism of Vandenbosch and Huizenga 


AA^/^afiE* - Bn) 


(^2^/an{E* - Bn) - 2^Jaf{E* - Bj) 


where a„ and af are the level density parameters at the equilibrium deformation and saddle 
point, respectively. Bn is the neutron separation energy, R/ is the hssion barrier and Kq = 
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/2mrQ. The fission barriers Bf are taken to be the snm of the Thomas-Fermi barrier [32] 
pins the shell correction term, 


Bf = Bf + Ushell (2) 

Usheii is taken to be the microscopic energy of the Finite Range Droplet Model (FRDM) [33] . 
Angnlar momentnm J is treated by redncing the available energy in Eq. [^by the rotational 
energy of the fissioning nnclens and scaling the Thomas-Fermi fission barrier with the 
Sierk barrier [33] ■ 

The fade-ont of the shell correction with increasing excitation energy is treated throngh 
the level density parameter following the method of Ignatynk et ah [35] . 

a(U) = a{l +f{U)5W/U) (3) 

where U is the excitation energy, 6W = M^xp^Z, A) A, a) is the difference between 

the experimental mass and the theoretical mass within the FRDM (the shell correction to 
the mass formula), and. 


f(U) = l-exp{-\U) (4) 

is a semi-empirical formula that drives the energy dependence of a. The asymptotic level 
density parameter a is given by. 


a = aA -|- /3A‘^^^s (5) 

where s is the surface on the nucleus in units of the equivalent-size sphere. The nuclear 
surface area S is estimated using the standard expansion of the nuclear radius in spheri¬ 
cal harmonics, which for symmetric deformations (as in a nucleus undergoing fission) and 
ignoring higher order terms, is given by. 


where 


S = AnRl 


2 . 

l + -a. 



02 — 



/92 


( 6 ) 

(7) 
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and (52 is the calculated quadrupole deformation of the nuclear ground state within the 
FRDM. We use the coefficients obtained with a realistic Wood-Saxon potential [33] . 


a = 0.073,/3 = 0.095,7 = 0.061 MeV^ 

The present simulations take the output of GRAZING in the form of the excitation energy 
E* distributions of primary products {Z^A) for each partial wave, which is converted into 
a discrete cumulative probability function, which in turn is used to numerically select an 
event with the generation of a single random number. Fig. [T] shows the simulated excitation 
energy distribution in the ^^®Xe+^°®Pb reaction for the partial wave leading to the highest 
cross section for producing primary product {N = 126) at E^,rn. = 423,450, 526 and 617 
MeV. The most probable E* is 6.0,8.4,19.2 and 30.0 MeV, for partial waves L = 64,156, 288 
and 394h, respectively. The average transferred angular momentum J at the most probable 
E* is 0.031,0.23,1.1 and 3.5h, respectively. 

For each initial event {Z,A,E*,J), F„/F/ is calculated using Eq. [^assuming a/ = an- The 
calculated F„/F f is tested with a random number to decide whether neutron evaporation 
or hssion happens. If hssion happens, the testing of the event is terminated. If neutron 
evaporation happens, A is decreased by one mass number and E* is reduced by {Bn + En), 
where Bn is the neutron binding energy and En is the neutron kinetic energy sampled 
randomly with a Maxwellian probability function of nuclear temperature T = a/ aE*, 

En = -T* [log(ri) + log(r 2 )] (8) 

where ri and r 2 are two independent random numbers. If J > 0, it is assumed the evaporated 
neutron carries Ih of angular momentum. The procedure is iterated until E* < Bn- 

Each simulation is performed with the standard set of parameters of GRAZING ng and 
the de-excitation part is simulated with 10^^ cascades in a High-Performance Gomputing 
Gluster using 40 nodes. This large number of cascades is necessary in order to simulate 
events with the lowest cross sections. The angular momentum transferred to the primary 
products is rather modest and it is therefore assumed that J = Oh in all simulations except 
where otherwise indicated. 

In what follows we refer to the simulations described in this section as GRAZING-F. 
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III. COMPARISON WITH EXPERIMENTAL DATA 


We have gathered an extensive set of experimental data to compare with simnlations. 
The reactions we have stndied can be divided into two categories; reactions in which the 
target is a Pb-like nnclens or is an actinide. In the former case, the hssility of the primary 
fragments is relatively low and hssion may be relevant only in the case of very high excitation 
energy. 

Some of the experimental stndies have been done with very thick targets, which pose 
a difhcnlty when comparing to simnlation since the reported cross section represents an 
integrated qnantity between the incident and exit projectile energies. If the projectile stops 
in the target material, the cross section represents an integrated qnantity down to the 
interaction barrier of the reaction. For thin-target experiments (for which the projectile exits 
the target), the projectile energy nsed in the simnlations was assnmed to be the effective 
mid-target projectile energy, estimated with range tables [36]. For thick-target experiments 
(for which the projectile stops in the target), the simnlations were done in snitable slices of 
the effective target thickness (the range np to the interaction barrier) and the cross section 
was calcnlated as the weighted mean of the slice cross section simnlated at the mid-slice 
energy. 

We have stndied only the yields of snrviving target-like prodncts. Table [T] lists the re¬ 
actions simnlated, the interaction barriers, the simnlated transfer and transfer-hssion cross 
sections. The last colnmn of the table is the reference to the experimental data nsed in the 
comparisons. 

A. The reactions 

The 238pj_|_238u 238 u_|_248qj^ reactions were stndied in the late 70s (Ref. [3|) and 1982 

(Ref. m), to determine the feasibility of nsing mnlti-nncleon transfer and deep-inelastic re¬ 
actions to synthesize snperheavy elements, ^sspj bombarded thick 238u 248Qjg^ 

targets and radiochemical methods were employed to dednce cross sections of actinide iso¬ 
topes. The experimental data was later reexamined by Kratz et al. [2Z] . The data reported 
in this latter paper form the basis of the present comparison with simnlations. The two 
systems have also been modeled within the diffnsion model |3H| and the dynamical model 
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based on the multi-dimensional Langevin equations |10j . 

The reaction was experimentally studied at entering projectile energy of 

1760 MeV with a target thick enough to stop the projectiles. The mid-target energy is 1650 
MeV. For the purpose of the simulations, the effective target thickness (4.8 mg/cm^) was 
divided in ten equal slices (equivalent to a stack of ten thin targets of 0.48 mg/cm^ each) 
and GRAZING was run at the equivalent mid-slice energy. The upper panel of Fig. shows 
a comparison between the cross section obtained by weighting the ten yields of Z = 98 
primary products (solid line) and the yield resulting from the effective mid-target energy 
alone (dashed line.) The weighted distribution is broader because it includes partial distri¬ 
butions at higher energies. The lower panel in Fig. shows the deviation. Assuming a single 
mid-target energy for this reaction may result in a systematical error of ~ 10% around the 
most probable mass, and more than 50% at the extremes. This result justihes the use of the 
weighted procedure at the expense of considerable computing time. 

Fig. i shows the direct comparison between experimental data and simulations with 
GRAZING-F for actinides with Z=97-101 in the ^^®U+^^®Gm reaction. The experimental 
data of Ref. [3] is plotted as solid symbols, the simulated yield of surviving products as 
solid lines and the primary product yields as dotted lines. The experimental cross section 
for ^^^Bk is a lower limit, which is denoted by upward arrow in the panel for Z = 97. The 
agreement for -I- Ip, +2p and -|-3p transfers is remarkable, whereas the simulation is less 
successful for larger p transfers. The reason for this discrepancy is that GRAZING predicts 
insufficient primary neutron transfers for these nuclei, as can be seen by comparing the 
simulated primary and secondary yields. The larger excitation energy moves the surviving 
product distribution towards lower A values. For the simulation to reproduce the yields of 
Fm and Md nuclei at least 2 and 4 additional neutrons would, on average, be required to be 
transferred to the primary products. The odd-even effects displayed by the simulations are 
a direct consequence of the odd-even effects introduced in the calculation of F„/F f in the 
de-excitation stage. The yields predicted by the Langevin model [lO] are shown as dashed 
lines in Fig. (see Fig. 4 in Ref. [TO].) 

The 238pj_j_238pj ];gaQtion was experimentally studied at four energies for which independent 
yields are reported. The targets used were thick enough to stop the projectiles. The entering 
energies were 2059 (11.6), 1785 (5.7), 1628 (2.4) and 1545 (0.65) MeV (mg/cm^), with mid¬ 
target energies of 1787, 1650, 1571 and 1530 MeV, respectively. Given in parenthesis is the 
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effective target thickness. For the purpose of simulations with GRAZING-F, the effective 
target thicknesses were divided in 20, 10, 5 and 1 shce(s), respectively. In the 2059 MeV 
reaction, the transfer-hssion cross section in the hrst slice is ~ 90% of the transfer cross 
section and decreases to ~ 4% in the last slice. Fig. shows the transfer and transfer-hssion 
cross section as a function of mid-slice projectile energy. The dependency of the transfer- 
hssion cross section on energy is determined by both the excitation energy distribution of 
the primary fragments and F„/F j. 

Fig. I shows the results for the 238pj_j_238pj r 0 action. The simulations do not reproduce 
the data. GRAZING does not predict > -|-5p transfers at the lowest energy (not shown 
in Fig. 1^ for that reason), and > +Qp transfers in the 1628 MeV reaction. As in the 
case of the reaction, GRAZING seems to underpredict the how of neutrons 

in > -|-4p transfers. In the 238pj_j_238pj taction, at least 5 additional neutrons would on 
average be required to be transferred to the primary products in order for GRAZING-F 
to reproduce the locations of the maximum yields. The yields predicted by the Langevin 
model HDl for the Eiab = 1785 MeV reaction are shown as dashed lines in Fig. Gomparing 
both simulations in the two reactions, we may conclude that the Langevin model seem to 
reproduce fairly well the yields of massive transfers (> +5p), as in the 238pj_j_238pj taction, 
whereas GRAZING-F seems to better reproduce the yields of a few-nucleon transfers (< 
-|-4p), as in the ^^®U-l-^^®Gm reaction (see Fig. [^) 

If we assume GRAZING-F is able to reproduce the yields of < +Ap transfers reasonably 
well, then GRAZING-F predicts substantial yields of unknown actinides in the 238pj_j_238pj 
reaction at Eiab = 2059 MeV. Fig. shows the predicted production cross sections oi Z = 
93 — 94 isotopes. Open circles represent unknown actinides. The predicted cross sections 
that are measurable (> 100 nb) are listed in Table [11} 

B. The reaction 

The ^^®Gm reactions were measured in order to study the influence of the 

projectile N/Z ratio in the production of actinides |H]. This work used a thin Gm target 
and the simulations were therefore done at the mid-target energy of the projectile. The 
cross sections in the ^^®Xe-t-^^®Gm reaction were measured by chemical separation with the 
intent to determine the formation cross sections of unknown actinides [6j. In this case the 


simulation was done at energy Eiab = 769, which we have estimated to be the effective 
mid-target energy for the reaction with entering energy of 790 MeV. 

Fig. 0 shows the comparison between experiment (solid symbols), the simulation of sec¬ 
ondary (solid lines) and primary product yields (dotted lines). The agreement between 
prediction and measurement is quite reasonable. GRAZING-F seems to do a fair job in 
predicting the cross sections of +2p transfer reactions but fails to reproduce the data for 
larger p transfers and some below-target yields. 

C. The ^^®Xe-|-^^®Cf reaction 

The ^^®Xe+^^®Gf reaction was measured with the intent to study the feasibility of us¬ 
ing low-energy multi-nucleon transfer reactions to produce new actinide and transactinide 
isotopes M- Fig. I shows a comparison between experimental data and simulation with 
GRAZING-F for the three mid-target energies studied, Eiab = 749, 813, 877 MeV, respec¬ 
tively. In this case, the simulations seem to predict the location of the maximum of the mass 
yields, but fail to predict the absolute values, overestimating the cross sections by an order 
of magnitude. 

D. The ^^®Xe-|-^^^Pu reaction 

The reaction ^^®Xe-l-^‘^'^Pu at Eiab = 835 MeV was used to produce and study the decay 
properties of the neutron-rich isotopes ^^^Np and ^"^^Np [7j. In Fig. [^we show the measured 
production cross section of Np isotopes compared to the predictions of GRAZING-F simula¬ 
tions. The simulations were done at energy Eiab = 826 MeV, which we have estimated to be 
the mid-target energy. The predicted yield pattern is more neutron rich than the observed 
yield pattern but is similar in shape. 

E. The ®®Kr-|-^^®Cm reaction 

The ®®Kr-|-^^®Gm reaction was studied experimentally in the 1980s [6]. Our simulations 
were done at Eiab = 435 and 457 MeV, corresponding to the entrance projectile energy, and 
Eiab = 520 MeV, which we have estimated to be the mid-target energy for the reaction with 
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entering energy of 546 MeV. The former two energies are either below or at the interaction 
barriers (see Table |^) Fig. 10 shows the comparison between experiment and simulations of 
secondary (solid lines) and primary product yields (dotted lines.) The observed yields are 
generally well represented by the GRAZING-F calculations. 


IV. PREDICTIONS 


A. The reaction 


2081 


The study of the ^^®Xe+^‘’®Pb reaction was hrst proposed by Zagrebaev and Greiner [39] 
as a way of demonstrating how nuclear structure effects could be influencing the flow of 
nucleons in low-energy multi-nucleon transfer reactions towards both the Z = 82 and N = 
126 closed shells. A dynamical model based on the multi-dimensional Langevin equations 
was used and the reaction was studied at Ec.m. = 450 MeV. Mass-energy distributions 
of the ^^®Xe+^°^Pb reaction have been measured recently with a double-arm time-of-flight 


spectrometer at Ec.m. = 423, 526 and 617 MeV M- In Fig. ^we show the yields from the 
GRAZING-F simulations at energies Ec.m. = 423,450,526 and 617 MeV for transfers where 
unknown N = 126 products are produced. (Unknown isotopes are plotted as open circles, 
while unknown N = 126 isotopes are plotted as solid circles.) The transfer-hssion cross 
section increases substantially in going from the lowest to the highest energy, from ~ 10 
mb to ~ 300 mb, as can be intuitively expected due to the larger excitation energy of the 
primary products. The simulated yields for yy^Ir and 7 g^Pt peak at Ec.m. ~ 750 MeV, with 
cross sections of 1.7 and 23 /rb, respectively. Assuming realistically a beam intensity of 100 
pnA and a target thickness of 1 mg/cm^, the production rates at this energy would be 3 
and 40 s“^, respectively. The range up to the interaction barrier is 2.75 mg/cm^. Hence, 
the maximum production rate at this energy and beam intensity would be 8 and 110 s“^, 
respectively, assuming the cross section between Unt and Ec.m. varies slowly with energy. 
The simulations suggest that the ^^®Xe -|- taction can be an important source of new 

neutron-rich nuclei near the N = 126 shell. 
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B. The reaction 


The reaction at Eiab = 9 MeV/A has been proposed as a = 126 “factory” 

based upon calculations using the GRAZING code without de-excitation by fission |TT]. If 
these predictions are correct, the properties of many unknown neutron-rich N = 126 nuclei 
could be studied with intense ^^®Xe beams and thick ^®®Pt targets. Although this reaction 
has not been studied experimentally, we have performed simulations with GRAZING-F 
in case fission plays a role in the de-excitation of primary reaction products. We find 
that GRAZING-F (assuming J = Oh) predicts a transfer-hssion cross section of ~ 30 mb. 
Gompared to the transfer cross section of ~ 5 b, fission does not seem to play a role if J is 
low. Even if the transferred angular momentum is large, say J = 30h, the largest angular 
momentum transferred predicted by GRAZING, the isotopic yields are essentially the same. 
Hence, fission competition nor angular momentum seem to play a role in this reaction. 


In Table III we show the maximum production rates for N = 126 isotopes by assuming 
a beam intensity of 1 p/iA and a target thickness equivalent to the range from the entrance 
energy to the interaction barrier. The simulations suggest that the use of ^®®Pt as & N = 126 
“factory” is justihed HU and may have a signihcant advantage over ^°®Pb, as the simulated 
transfer cross section in the former case may be a factor of two higher. 


C. The ^^^Xe-|-^^®Cm reaction 


The ^^^Xe-l-^^®Gm reaction at Eiab = 800 MeV has been proposed as an example reaction 
to be studied at EURISOL in a series of meetings and workshops [l2lll3]. The presentations 
suggest that GRAZING predicts this reaction has production cross sections of the order of 
0.1 mb for U, 1 mb for Np and Pu, and 10 mb for Am neutron-rich isotopes. However, these 
calculations did not consider neutron decay and thus represent yields of primary fragments 


only H- Fig. [l^ shows the predictions of GRAZING-F. Unknown isotopes are shown as 
open circles. The present simulations predict cross sections of the order of /ib to nb for the 
most neutron-rich unknown actinides. The transfer cross section is estimated to be ~ 7 b, 
whereas the transfer-fission cross section is estimated to be ~ 160 mb indicating that fission 
is not an important decay mode for the most n-rich products. GRAZING does not predict 
the production of Z = 92 isotopes. 


11 




The trend as a function of projectile mass is shown in Fig. 1^, where we plot the cross sec¬ 
tion of surviving nuclei of —2p, —Ip, Op, -|-lp and +2p transfers for the i29,i32,i36,i44^g_|_248Qj^ 
at Ec,m./Vc = 1-05. Unknown isotopes are shown as open symbols. If we focus on cases 
that are well described by GRAZING-F {+2p transfers), these simulations seem to indicate 
that a more neutron-rich projectile does indeed produce more neutron-rich products, but 
the advantage of going from the most neutron-rich stable Xe to a neutron-rich radioactive 
Xe projectile may not be as pronounced as claimed or hoped. One notes furthermore that 
the projected intensities of ^^^Xe beams at modern radioactive beam facilities are a tiny 
fraction of the intensities of the stable Xe beams. 


D. The ^^Kr-l-2^^Cm reaction 


The ®"^Kr-|-^^®Gm reaction has been simulated at Ef.,m./Vc = 1-45 and compared to the 
®®Kr-|-^^®Gm reaction in Fig. 14 Unknown isotopes are shown as open circles. The ®^Kr 


simulations predict substantial cross sections for unknown neutron-rich actinide isotopes 
compared to ®®Kr. For example, the predicted production cross section for ^“^^Pu is ~ 0.5 
mb in the ®^Kr induced reaction, whereas the cross section is ~ 0.02 mb in the ®®Kr induced 
reaction. Nonetheless, current ®"^Kr intensities are far too low for this reaction to be feasible 
to produce actinides. For example, the maximum production rate of ^"^^Pu would be ~ 15 
per year assuming the current GARIBU beam intensity estimate of 15 s“^. 

The simulations also predict the absence of larger p transfers (> -l-5p) in the ^'^Kr induced 
reaction. 


E. The reaction 

The reaction has been suggested to be studied with a mass separator like the 

Fragment Mass Analyzer (FMA) at Argonne National Laboratory. The reason is that there 
is some evidence that large yields of transfer products could be observed close to 0° |13]. If 
this is the case, the yields of short-lived neutron-rich actinides could be measured and the 
theory of Zagrebaev and Greiner [10] could be tested. The particular choice, ^^®U-|-^^®Bk, 
has been suggested because of convenience; ^^®U is the heaviest projectile accelerated by 
ATLAS and a thin ^^®Bk target has recently become available. 
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Fig. shows the cross section of surviving nuclei predicted by GRAZING-F when 
Ec.m./yc = 1-05 and 1.45, respectively. Unknown isotopes are shown as open circles. The 
simulations predict substantial cross sections for unknown U and Pu isotopes at both en¬ 
ergies, with the larger cross sections associated with the low-energy reaction due to hssion 
competition. Table |IV| shows the yields of unknown U and Pu isotopes assuming a beam 
intensity of 100 pnA, a target thickness of 0.3 mg/cm^ and one day of irradiation. In the 
high-energy reaction GRAZING-F predicts the production of unknown neutron-rich Md, No 
and Lr isotopes with cross sections below 1 nb. Under the above assumptions, ^®^Md would 
have a yield of ~ 30, ^®^No ~ 10 and ^®^Lr ~ 1 nuclei, respectively. 

V. CONCLUSIONS AND DISCUSSION 

From the comparison of available experimental data with simulations of GRAZING-F we 
conclude that it is able to reproduce the yields of above-target products of +lp, +2p and 
-|-3p transfers reasonably well. The yields of -|-lp, +2p and -|-3p transfers in the ^^®U-f^^®Gm 
reaction, for example, are exceptionally well reproduced. The yields of products involving 
larger proton transfers (> -|-3p) start to deviate substantially from experimental data pri¬ 
marily because GRAZING seems to underestimate the flow of neutrons. The usefulness 
of very neutron-rich radioactive beams, however intense, is predicted to be doubtful com¬ 
pared to neutron-rich stable beams, e.g. ®®Kr and ^^®Xe, which are far more intense than 
any predicted intensities of the radioactive beams at planned radioactive beam facilities. 
GRAZING-F predicts substantial yields of unknown actinide isotopes under special condi¬ 
tions. For example, in the low energy U-|-Bk reaction, several unknown U and Pu isotopes 
could be produced with measurable yields. The production of unknown N = 126 isotopes 
is predicted to be better accomplished by the use of the Xe+Pt reaction. The low transfer 
cross section and the high transfer-hssion cross section associated with the Xe-|-Pb reaction 
makes it a less attractive candidate. Of course, these are predictions that must be verihed 
by experiment. 
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FIG. 1: Excitation energy distributions predicted by GRAZING in the reaction ^^®Xe+^®®Pb at 
Ec.m. = 423,450, 526 and 617 MeV for the partial wave leading to the highest cross section for 
producing primary product 78 ‘^Pti 26 - The partial wave L is given in the panels. 
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TABLE I: List of reactions simulated with GRAZING-F. The Coulomb and interaction barriers 
are calculated within the Bass model [l6]. Cross sections are taken from simulations. The last 
column gives the reference to the experimental data if it exists. 


Reaction Eiab (MeV) E,.m Vc (MeV) Vint (MeV) Ec.m/Vc (mb) (“b) Ref. 


136xe+208pb 

701.4 

423.0 

423.5 

430.5 

1.00 

2010 

7 

m 


746.3 

450.0 

423.5 

430.5 

1.06 

2340 

29 



872.8 

526.0 

423.5 

430.5 

1.24 

2700 

122 

m 


1024.3 

617.0 

423.5 

430.5 

1.46 

2900 

268 

m 

136xe+198pt 

1224.0 

604.9 

405.9 

412.1 

1.49 

5340 

29 


86Kr+248cm 

435.0 

323.0 

340.4 

344.3 

0.95 

6360 

46 

m 


457.0 

339.3 

340.4 

344.3 

1.00 

6590 

210 

m 


520.0 

385.1 

340.4 

344.3 

1.13 

7000 

3880 

m 


667.4 

494.0 

340.4 

344.3 

1.45 

7540 

6700 


94Kr+248Cj^ 

677.8 

490.0 

336.9 

340.7 

1.45 

8300 

7450 


136Xe+244pu 

826.0 

528.8 

473.9 

482.7 

1.12 

2750 

1510 

m 

^^®Xe+^^®Cm 

780.0 

511.6 

486.0 

495.3 

1.05 

7330 

860 

m 

^^^Xe+^^®Cm 

782.0 

508.9 

484.6 

493.8 

1.05 

7090 

670 



805.0 

523.8 

484.6 

493.8 

1.08 

7300 

800 

0 

136xe+248Cm 

769.0 

496.6 

482.7 

492.0 

1.03 

7330 

260 

m 


785.0 

505.6 

482.7 

492.0 

1.05 

7050 

810 


144Xe+248Cm 

800.0 

504.7 

479.2 

488.3 

1.05 

7150 

160 


136Xe+249cf 

749.0 

483.1 

492.5 

502.1 

0.98 

2290 

1680 

0 


813.0 

524.3 

492.5 

502.1 

1.06 

2710 

1430 

m 


877.0 

565.4 

492.5 

502.1 

1.15 

2910 

1670 

0 

238p_j_238p 

2059.0 

1024.8 

735.2 

754.6 

1.39 

9310 

4710 

m 


1785.0 

899.9 

735.2 

754.6 

1.21 

9060 

1760 

m 


1628.0 

811.0 

735.2 

754.6 

1.10 

8930 

750 

m 


1545.0 

769.8 

735.2 

754.6 

1.05 

8750 

350 

m 

238u+248Cm 

1760.0 

894.6 

762.6 

783.2 

1.17 

8860 

1400 

m 

238u+249Bk 

1587.9 

809.0 

770.1 

791.0 

1.05 

8780 

6500 



2195.4 

1117.0 

770.1 

791.0 

1.45 

9980 

7960 
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TABLE II: Predicted cross sections of unknown U and Np isotopes in the 238-y_j_238pj 


Eiab = 2059 MeV. Only cross sections > 100 pb are listed. 


Isotope 

fj ifih) 

243u 

21.8 

244u 

50.8 

245u 

3.7 

246u 

2.4 

247u 

0.25 

245Np 

2.9 

246 p^p 

0.22 

247Np 

0.35 


TABLE III: Maximum production rates of Ai = 126 isotopes in the ^^®Xe+^®®Pt reaction at Eiab = 9 
MeV/A simulated with GRAZING-F assuming a beam current of 1 P/uA and a target thickness 
equivalent to the range from the entrance energy to the interaction barrier. 


Isotope Rmax (s 


2.6 X 10^ 


4.7 X 10® 

202q 

76 

5.5 X 10^ 

201r„ 

75 

4.0 X 10^ 


20 



TABLE IV: Yields V of U and Pu isotopes in the reaction at Ec.m. = 809 MeV 

simulated with GRAZING-F assuming a beam current of 100 pnA, a target thickness of 0.3 mg/cm^ 
and 1 day irradiation. 


Isotope 

N 

Isotope 


N 

244u 

1.9 

X 

10^ 

248pu 

1.6 

X 

10^ 

245u 

5.7 

X 

10^ 

249pu 

2.8 

X 

10® 

246u 

1.3 

X 

10^ 

250pu 

3.2 

X 

10® 

247u 

2.9 

X 

10^ 

251pu 

to 

o 

X 

10® 

248u 

2.8 

X 

10^ 

252pu 

CO 

o 

X 

10® 

249u 

6.7 

X 

10^ 

253pu 

7.0 

X 

10^ 

250^ 

2.5 

X 

10^ 

254pu 

1.4 

X 

10^ 
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FIG. 2; Panel a) shows a comparison between simulations of the primary yields of Z = 98 products 
in the reaction with GRAZING by constructing a weighted cross-section (solid line) 

and a single simulation at the mid-target energy (dashed line.) The effective target thickness is 
4.8 mg/cm^ and the weighted simulation is made by assuming a stack of ten identical target slices. 
Panel b) shows the deviation. 
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FIG. 3: Cross sections of surviving nuclei in the reaction at entering projectile energy 

Eiab = 1760 MeV. Experimental data from Ref. m are shown as solid symbols and the predicted 
cross sections (GRAZING-F) as solid lines. The measurement of is a lower limit which is 

indicated with an arrow. Dotted lines show the predicted yields of primary products. Dashed lines 
shows the Langevin simulations in Ref. uni 
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FIG. 4: Predicted transfer (GRAZING) and transfer-fission cross sections (GRAZING-F) as a 
function of mid-slice energy in the 238-y_|_238'[j i-gag^ion at entering energy Eiab = 2059 MeV. 
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FIG. 5: Cross sections of surviving nuclei in the reaction 238-jj_j_238u _ 3 ^ 028 ,1785 and 2059 

MeV. Experimental data from Ref. m are shown as solid symbols and predicted cross sections 
(GRAZING-F) as solid lines. Dotted lines show the predicted yield of primary products. Dashed 
lines shows the Langevin simulations in Ref. [lOj . 
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FIG. 6: Predicted (GRAZING-F) cross sections of surviving nuclei with Z = 92 — 94 in the 
238 u_|_238u reaction at entering energy Eiab = 2059 MeV. Unknown isotopes are shown as open 
circles. 
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FIG. 7: Predicted (GRAZING-F) cross sections of surviving nuclei in the reaction ^^®Xe+^^®Gm 
at Eiab = 780 MeV, at Eiab = 805 MeV and at Eiab = 769 MeV (solid 

lines) compared to experimental data [HlIH] (solid symbols). Predicted primary product yields are 
shown as dotted lines. 
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FIG. 8; Predicted (GRAZING-F) cross sections of surviving nuclei in the reaction ^^®Xe+^^®Gf at 
Eiab = 749,813,877 MeV (solid lines) compared to experimental data [9] (solid symbols). 










A 


FIG. 9: Predicted (GRAZING-F) cross sections of surviving nuclei in the reaction ^^^Xe+^'^^Pu at 
Eiab = 826 MeV (solid lines) compared to experimental data [7] (solid symbols). Predicted primary 
product yields are shown as dotted lines. 
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FIG. 10: Predicted (GRAZING-F) cross sections of surviving nuclei in the reaction ®®Kr+^^®Gm 
at Eiab = 435,457 and 520 MeV (solid lines) compared to experimental data [6] (solid symbols). 
Predicted primary product yields are shown as dotted lines. 
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FIG. 11; Cross sections of surviving nuclei in the reaction ^^®Xe+^®®Pb at Ec.m. = 423,450,526 
and 617 MeV. The simulations (GRAZING-F) are shown as solid lines. Unknown isotopes are 
shown as open circles and unknown isotopes with N = 126 are shown as solid circles. 
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FIG. 13: Cross sections of surviving nuclei in the reaction 129,132,i36,i44^g_j_248Qj^ = 

1.05. The predictions (GRAZING-F) are shown as solid lines. Unknown isotopes are shown as 
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open circles. 
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